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The 71 ^ — rj — 1]' mixing in a generalized multi-quark interaction scheme. 
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We investigate the isospin symmetry breaking effects within a recently derived Nambu-Jona- 
Lasinio related model by fitting the measured pseudoscalar meson masses and weak decay couplings 
/tt, /x. Our model contains the next to leading order terms in the 1/Nc expansion of the effective 
multi-quark Lagrangian, including the ones that break the chiral symmetry explicitly. We show the 
important phenomenological role of these interactions: (1) they lead to an accurate fit of the low- 
lying pseudoscalar nonet characteristics; (2) they account for a very good agreement of the current 
quark masses with the present PDG values; (3) they reduce by 40% the ratio e/e' of the — rj 
and 7r° — rj' mixing angles, as compared to the case that contemplates explicit breaking only in the 
leading order, bringing it in consonance with the quoted values in the literature. The conventional 
NJL-type models fail in the joint description of these parameters. 

PACS numbers: 11.30.Rd, ll.30.Qc, 12.39.Fe, 12.40.Yx, 14.40.Aq, 14.65.Bt 


I. Introduction 

The strong isospin symmetry is considered to be a 
very good approximation in the empirical description of 
a large bulk of strong interaction processes. This is re¬ 
lated to the hierarchy in which breaking of the chiral 
symmetry SU{3)l x SU{ 3 )ii by different current quark 
masses occurs, down to SU{2)i x U{1)y flavor symme¬ 
try if mu,rnd mg. In the case of the pseudoscalar 
mesons it is accurate at the order of the ratio of the 
light and strange current quark masses (m„ — md)/ms 
[lj,[3 and explains partly the small meson mass differ¬ 
ences within charged isospin multiplets. A further source 
of isospin breaking is due to the electromagnetic interac¬ 
tions, which are expected to be suppressed at the scale 
of strong interactions. 

A detailed quantitative analysis however requires 
isospin breaking corrections to be taken into account in a 
series of low energy phenomena, such as: the description 
of mass splittings of mesons, and Dashen’s theorem 3; 
sum rules for quark condensates 0,0; kaon decays ^ i 
TT — TT i,0 and TT — K scattering 0,0 in relation to 
mesonic atoms, p and T decays involving 77(77') 

mesons [l ^ , ao — fo mixing [l^ in the scalar meson sec¬ 
tor. 

Strong isospin breaking effects become particularly rel¬ 
evant if a certain process depends crucially on the differ¬ 
ences of the light quark masses. If in addition the elec¬ 
tromagnetic interactions are a subleading effect, these 
processes provide for ideal tools in a quantitative analy¬ 
sis of quark mass ratios. In the latter category are the 
77, 77' — >■ Stt decays, the 7r° — 77 and tt^ — 77' mixings, as well 
as the p — oj mixing in the vector channels. 

Isospin breaking associated with the 7r° — 77 — 77' sys¬ 
tem has long been known to play a role in the Standard 
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Model prediction of the CP violation related ratio {t)cp 
[mil representing a substantial correction to the QCD 
penguin contributions [l^. It affects the value of the 
^ 7r°7r° transition through the dominant QCD Qq 
penguin operator, which is one of the sources of uncer¬ 
tainties in the determination of {^)cp ^ recent 

review see [10 . 

In chiral perturbation theory (ChPT)[ 23 ,[iil ,0 the 
7r° — 77 mixing angle occurs already at order and was 
first evaluated to order in the context of Ki^ form 
factors in pTI ]. 

In the analysis of 77 — 77' mixing of [ 2 ^ the U{1)a 
anomaly is described by the gluon transition matrix el¬ 
ement < 0|^GG|77i > and the quark flavor basis has 
been used. As shown in [ 2 ^ this basis is favored, as one 
of the two mixing angles is indeed small. The decay con¬ 
stants follow the pattern of particle state mixing in that 
basis. It has been shown that this approach leads to re¬ 
sults consistent with many observables related to 77 — 77' 
mixing In it has been extended to include 

the mixing to the neutral pion. 

In the present work we address the tt^ — rj — rj' mixings 
resulting from a recently proposed Lagrangian 
which is reviewed below. In this effective Lagrangian 
approach built from all spin 0 and non-derivative multi¬ 
quark interactions relevant at the scale of spontaneous 
chiral symmetry breaking, the complete set of interac¬ 
tions which break explicitely the chiral symmetry was 
included for the first time. This Lagrangian repre¬ 
sents a generalization of the original Nambu-Jona-Lasinio 
[H, [ 2 ^ model extended to the realistic three flavor and 
color case with U{ 1 )a breaking six-quark ’t Hooft in¬ 
teractions l30l-l^ and an appropriate set of eight-quark 
interactions [45l| . The last ones complete the number 
of vertices which are important in four dimensions for 
dynamical SU{ 3 )l x SU{ 3 )ji chiral symmetry breaking 
[4^ 1^ . The Lagrangian considers all interactions rele¬ 
vant at the same order in large Nc counting as the U{ 1 )a 
anomaly term. 

The role of the new interactions contained in the ex¬ 
plicit symmetry breaking (ESB) vertices has been an- 
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alyzed at meson tree level approximation and in the 
isospin limit in connection with the low lying characteris¬ 
tics of the pseudoscalar and scalar meson nonets , [l^ 
and in the T — fj, phase diagram associated with chiral 
transitions [i^. An unprecedented accuracy for the de¬ 
scription of the spectra has been achieved. One should 
stress that the present Lagrangian is able to account 
properly for the SU{3) breaking effects in the descrip¬ 
tion of the weak decay constants /tt and /k, in addition 
to yield the correct empirical 77 , 77 ' and K meson masses, 
as well as the anomalous two photon decays of 71,77,77', 
in an unified description, which was an open problem for 
model versions without the ESB terms. 

The paper is organized as follows. In the next section 
is presented the effective multiquark Lagrangian and its 
bosonized form, the associated Nc counting is reviewed, 
in section III we address the mixing in the tt — r] — rj' 
system, the choices of representation of states, the decay 
parameters in the flavor basis, and the compliance of the 
model in the approximation considered with the decay 
parameters transforming as the states. In section IV we 
present and discuss the numerical fits of the mass spectra 
and decay parameters. We end with a summary of the 
main results. 

II. Survey of the model Lagrangian 

II a. Multiquark picture 

The Langrangian considered is built from all spin 0 
chiral SU{3)l x SU{3)}i symmetric and parity conserv¬ 
ing combinations relevant at the scale A of spontaneous 
breaking of chiral symmetry. This means that in the cor¬ 
responding effective potential are kept all the interactions 
which do not vanish in the A —^ 00 limit. These con¬ 
sist of vertices involving non-derivative quark-antiquark 
fields, denoted by E, with E = ^{sa — ipa)^a and 
Sa = Pa = 9*75Aa9, Aa being the sta ndard Gell- 

Mann matrices for a = 1...8 and Aq = -^2/3 x 1, as well 
as of their interactions with the external sources x which 
transform as y = (3,3*) under SU{3)l x SU{3)h, 

L = qix'"df,q +L,nt +Lx- ( 1 ) 

The term 

L^nt = ^tr(EtE) + ^(detE + detEt) 

+ ||(trEtE)%||tr(EtEEtE). (2) 

is well known. Here and elsewhere we use barred quan¬ 
tities for any dimensionless coupling, these are related 
to the dimensionful ones through powers of the scale A, 
Pi = ffi/A'’'. Lint contains the leading order (LO) in Nc 
four quark (q) NJL interactions with coupling G, gener¬ 
alized to the 3 flavor case, the NLO 6q ’t Hooft U{ 1 )a 
breaking flavor determinant with coupling k, and the two 
possible 8q interactions ~ 51,521 which have the same Nc 
counting as the ’t Hooft term. We refer to [2^,[23 for 


a detailed discussion of the large Nc counting scheme 
which complies with the counting rules based on powers 
of the scale A of spontaneous breaking of chiral symme¬ 
try. Rephrasing it, means that the terms which vanish 
as A —>■ 00 are the same that also vanish on grounds of 
Nc — >■ 00. Applying these rules, summarized after eq. SI 
the following set of terms involving interactions with the 
sources y emerge, which act also up to the same order in 
Nc counting as the’t Hooft term 
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Lx = ^ Li, (3) 

i=0 

where 

Lo = -tr(E’''y-ty'''E) 

^1 — “t“ h.C. 

^2 — 'j^Gijk^mnlXi'rn^jn^kl h.C. 

L 3 = ||tr (EtEEty) + h.C. 

Li = ^tr (E’I'E) tr (E'I'y)-h h.C. 

Ls = ||tr (El'yEty) -p h.c. 

Le = ||tr (EEtyyt + 

Lr = (trE'l'y-h/i.c.)^ 

Ls = || (trE'I'y-h.c.)^ 

Lg = -^tr(E'l'yy'l'y)-p/i.c. 

Lio = -^tr (y'I'y) tr (y'I'E)-h h.c. (4) 

We recall that the Nc book keeping is as follows: 

E~V„ A~fV0 = l, G^l/Nc, n^l/Nl 
ki,59.5io ^ 1/A^c, 

K 2 , 35,56,57,58 ^ 53,54 ^ 1/A^c- (5) 

The scale A which gives the right dimensionality to 
the multiquark vertices is related with the cutoff on the 
divergent quark one-loop integrals, its Nc counting is dic¬ 
tated by the gap equations in the phase of spontaneously 
broken chiral symmetry, eqs. (I2l]),(l22l),(l23l) below. 

The terms L 1 ...L 1 Q of LaKrangi an ffl have been intro¬ 
duced recently in the model [^,[^ and generalize to 
NLO in Nc the explicit symmetry breaking (ESB) stan¬ 
dard LO mass term Lq. 

Since all the blocks Lq...Liq conform with the sym¬ 
metry pattern of the model one is free to choose for the 
source the constant valued matrix y = A4/2, 

M = /iaAa = diag{pLu,Pd,Ps)- 


( 6 ) 
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TABLE I: The pseudoscalar masses and weak decay constans (all in MeV) in the isospin limit used as input (marked with *) 
for different sets of the model. Parameter sets (a),(b) contain explicit symmetry breaking interactions (see Table HUll a nd allow 
for a fit of the scalar masses and strong decays as well, rria = 550 MeV, = 850 MeV, ma^ = = 980 MeV [^; set (c) 

does not. Set (a) corresponds to an octet-singlet mixing angle in the scalar sector of 6s = 27.5°, set (b) to 6s = 25°. 


Sets m.n- 

777 Tf 

nir, 


rrin' 

U 

/k 

a 138* 

494* 

547* 


958* 

92* 

113* 

b 138* 

494* 

547* 


958* 

92* 

113* 

c 138* 

494* 

475* 


958* 

92* 

115.7 

TABLE II: Parameter sets of the model: = 

= rud = m,ms, and A are given in MeV. The couplings have the following units: 

[G] = GeV-2, [ k ] = GeV®, [g^ 

l] ~ [52] = GeV ®. We also show here the values of constituent quark masses Mu = Md = M 

and Mb in MeV. See also caption of Table [H 






Sets m iria 

M 

Mb 

A 

G 

— K 

51 92 

a 4.0* 100* 

373 

544 

828 

10.48 

122 . 

3284 173 

b 4.0* 100* 

372 

542 

829 

9.83 

118.5 

3305 -158 

c 6.1 190 

375 

569 

836 

9.79 

138.2 

2500* 100* 

TABLE III: Explicit symmetry breaking interaction couplings. 

The couplings have the following units: [772] = GeV bs] = 

[qa] = GeV [55] = be] = [97] 

= bs] = GeV- 

See also caption of Tabled 



Sets k2 

-93 

54 

55 

-96 

-97 

98 

a 6.17 

6497 

1235 

213 

1642 

13.3 

-64 

b 5.61 

6472 

702 

210 

1668 

100 

-38 

c 0 

0 

0 

0 

0 

0 

0 

TABLE IV: The mixing angles 

in the 77 — 77' 

system in isospin limit, and related weak decay constants for the sets discussed 

and in comparison with different approaches. 

see also main text (the 

systematic error estimates given in 

[ 81(,[83 have been 

omitted here ). 







Sets 


do 



io. 

U 

is. 

a 

- 12 * 

-1.42 


-21.37 

1.172 

1.318 

b 

-15* 

-4.42 


-24.37 

1.172 

1.322 

c 

-14.5 

-2.82 


-24.78 

1.197 

1.365 

[ 2 ^ phen. 

-13.3 

- 6.8 


-19.4 

1.10 

1.19 

[ 2 ^ phen. 

-15.4 

-9.2 


- 21.2 

1.17 

1.26 

[7^ GHPT 

-10.5 

-1.5 


- 20.0 

1.24 

1.31 

GHPT 

- 

-4. 


-20.5 

1.10 

1.28 

sum rules 

- 

-15.6 


- 10.8 

1.39 

1.39 

Pade approximants ( 77 ) 

-16.4 

-11.3 


-21.3 

1.15 

1.22 

Pade approximants (77') 

-13.3 

-1.5 


-24.2 

1.28 

1.46 

BABAR (from ( 77 )) 

-21.7 

-26.7 


-16.5 

1.04 

0.98 

[8^ BABAR (from [ 8 IJ ( 77 ')) 

-17.7 

-15.6 


-19.9 

1.14 

1.11 

TABLE V: Empirical input used in the fits with isospin breaking, sets A and B with ESB interactions, set C without. Primes 

indicate which masses of the pion and kaon multiplets have been used for the fit, the other being output. 

Masses in units of 

MeV, angle tp in degrees. 







Sets 777° TTlJ 

rrin 

m'^ 

m\ 

777 ± 

U Jk 


A,B 136’ 136.6 

547 

958 

500 

494’ 

92 113 

39.7 

C 136’ 137.0 

477 

958 

501 

497’ 

92 116 

39.7 
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TABLE VI: Parameter sets with isospin breaking, set A with ESB interactions, set B without. The couplings have the following 
units: [G] = GeV^, [k] = GeV®, [ 51 ] = [ 52 ] = GeV®, [kj] = GeV®, [ 53 ] = [flr] = GeV®, [ 35 ] = [fle] = [flr] = [fls] = GeV^. 
A is given in MeV. See also caption of Table H] 


Sets 

G 

— K 

9i 

92 

A 

1^2 

93 

94 

95 

56 

97 

58 

A 

10.48 

116.8 

3284 

1237 

828.5 

6.24 

2365 

1182 

160 

712 

580 

44 

B 

10.48 

116.8 

3284 

1252 

828.5 

6.26 

2481 

1182 

151 

745 

591 

49 

G 

9.79 

137.4 

2500 

117 

835.7 

0 

0 

0 

0 

0 

0 

0 


TABLE 

and TT® 

1 VII: Isospin breaking parameter r = 
— rj, Tv^ — rj' mixing angles e and e . 

V+rV ’ cnrrent and constituent quark 

masses rriu^ 

md,ms Mu, Md, 

Ms in MeV 

Sets 

r 

mu 

md 

ms Mu 

Md 

Ms 

€ 

t 

€ 

€ 

A 

0.372* 

2.179 

4.760 

95* 372 

375 

544 

0.014* 

0.0037* 

3.78 

B 

0.372* 

2.166 

4.733 

95* 372 

375 

544 

0.017* 

0.0045* 

3.95 

G 

0.372* 

3.774 

8.246 

194 373 

380 

573 

0.022 

0.0025 

8.78 


Whenever convenient we use in the following the equiv¬ 
alent redefinition of flavor indices from a = 0,3,8 to 
i = u,d,s for any observable A [i^ 


A 


a 


5 



^ y/2 y/2 y/2\ 

Vi -Vi 0 

V 1 1 - 2 ] 


(7) 


The terms £i,£ R,£in are related to the Kaplan- 
Manohar ambiguity (SOl - l^ in the definition of the quark 
mass. In the present case this corresponds to the follow¬ 
ing freedom in transforming the external source 


X 


(c.) 


= X + ^ (det x't') X (x'^'x) ^ + ^XX^X 


with three independent constants Ci , which has the same 
symmetry transformation property as x- The description 
in terms of x^°‘^ in the place of x is equivalent, leading 
to the same Lagrangian, with some of the couplings re¬ 
defined as 

_ _ Cl _ _/ _ _ 

Ki= Kl + y, 95 ^ 95 = 95 - K 2 C 1 , 

- V , ^2 - 

97 ^ 97 — 97 + 98 ^ 9s — 98 + 

99^ 9g = 99 + C2- 2kiCi, 

510 g'lQ = gio -b C 3 -b 2 kiCi. (9) 

Note that the redefinition of the x fields ([5]) in the model 
contains more terms than the one usually considered, in 
the context of second order ChPT, which leads to the 
following redefinition of the current quark mass [50j | 

M^M{X)=M+XM{M^M)-^detM^. (10) 


It has been reported that the Kaplan-Manohar ambi¬ 
guity may be in conflict with the large N^. counting rules 
of ChPT [ 53 . Following from a threefold chiral expan¬ 
sion in the number of derivatives, powers of quark masses 


and powers of l/Nc, as well as the counting associated 
with the 9 angle related with the Cyi(l) sector, A is found 
to be suppressed to all orders in the large W limit [53 |. 

Regarding our model Lagrangian, the new couplings, 
the primed ones in (jOj, must not dominate over the 
dependence of the unprimed ones. This means that in 
their leading contributions they scale at most as the un¬ 
primed couplings, Ci ^ l/Nc, see JS]) . Attributing this 
counting to the Ci suffices to warrant that one can use 
the reparametrization freedom 0 , in particular to ob¬ 
tain k[ = g'g = Vio = 0- Indeed, remembering that our 
model Lagrangian has been constructed to incorporate 
the classes of multiquark interactions which do not vanish 
in the effective potential as Nc —>■ 00 in the phase of spon¬ 
taneously broken chiral symmetry, one gets for this case 
that Ki = — ci/ 2 , and that up to subleading corrections 
(not considered at this order of the effective potential), 
99 = -C2,510 = -C 3 , §5 = 95,97 = 97,98 = 98- 


II b. Bosonized version 


The low energy meson characteristica are obtained af¬ 
ter path integral bosonization of the quark Lagrangian 
(ED- Following one may equivalently use the intro¬ 
duced Sa,Pa as auxiliary fields, and a further set of phys¬ 
ical scalar and pseudoscalar fields a = GaXa, (j) = (paXa to 
obtain the vacuum persistence amplitude of the theory 
as 


Z= VqVqll GXp (i I d xLq (^, (7, 0) 


+ 00 

X / n n 



( 11 ) 


In these variables the Lagrangian reads 
L = Lq -\~ Laux 
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TABLE VIII: e and e' values in the literature. 



e 

/ 

e 

e 

V 

phen. 

0.014 

0.0037 

3.78 

[25] phen. 

0.017 ± 0.002 

0.004 ± 0.001 

4.25 ± 1.17 

[7^ ChPt NLO 

0.014^0.016 

- 

- 

[84] phen. 

0.021 

- 

- 

Exp. 

0.030 ± 0.002 

- 

- 

[86] Exp. 

0.026 ± 0.007 

- 

- 


Lq = q - (cr + M) - 775(/») q 

^aux — 5a(*^a “1“ ^a) “t“ Pa^a p) 

8 

+ '^ms,p,m). (12) 

i=2 


Here L' contains the ESB terms and Lq appears as SaVia 
in Laux- The external scalar fields a have been shifted to 
cr —7> cr + M, so that the expectation valne of the shifted 
fields in the vacuum corresponding to dynamically bro¬ 
ken chiral symmetry vanish. The expectation value of the 
unshifted scalar field < tr >= MaXa = diag{Mu, Md, Mg) 
corresponds to the point where the effective potential of 
the theory achieves its minimum, with M being the con¬ 
stituent quark masses. In (1121) the bilinear form in the 
quark fields Lq can be integrated out from the path inte¬ 
gral, and results in the fermion determinant (see m 
below), which generates the kinetic terms of the a,(j) 
fields. The remaining integrations are over a static non- 
Gaussian system of s,p fields, and are done in the sta¬ 
tionary phase approximation (SPA). First we present the 
results for the SPA integration, obtained from the ex¬ 
tremum conditions 


dL 

dSa 


= 0 , 


dL 

dpa 


= 0 , 


(13) 


which must be fulfilled in the neighbourhood of the uni¬ 
form vacuum state of the theory. The solutions of eq. 
m are seeked in the form: 


sf = ha + ab + h^abc^bCTc + h^abl(l^b(l3c + ■ ■ ■ 

Pa = h^lUb + h^abi<l^b<Jc + ■ ■ ■ (14) 


Eqs. m determine all coefficients of this expansion giv¬ 
ing rise to a system of cubic equations to obtain ha, eq. 
(di), and a full set of recurrence relations to find higher 
order coefficients in (ITdll . The result is cast in the form 


Laux = haaa + ^ h^ab^a^Tb + ^ h^ab 4'a(l 


(15) 


1 


^abl^bCTc 


' ( 2 ) 
^abc 


i(3) 

^bca 


^ (I^b4^c 


Here ha are related to the quark condensates, h^^J, 
contribute to the masses of scalar and pseudoscalar states 
respectively, and higher indexed h’s are the couplings 


that measure the strength of the meson-meson interac¬ 
tions. From eq. (USD and using © one obtains the fol¬ 
lowing system of cubic equations for the one index coef¬ 
ficients hi, {i = {u, d, s}) 

Mi — rrii + —tijkhjhk + — (2G + gih^ + gimh) + ^h\ 

+ ^ [^gshi -b gAh^ + 2(^5 -b ge)mihi + Agpmh] 
+K2tijkmjhk = 0. (16) 

Here Ujk is a totally symmetric quantity, whose nonzero 
components are tuds = 1; there is no summation over 
the open index i but we sum over the dummy indices, 
e.g. h"^ = hl^ + h"^ + h1,mh = ruahu + mdhd + rnghg. 
Regarding the h coefficients with more than one index, 
we indicate explicitly only the expression needed for the 
present study of the pseudoscalar masses (the complete 
expressions up to 3 indices can be found in 0,123) 

“2 {h^ab') = (2^ + gih^ + gAmh) Sab 

^-^abc {f^hc ‘2K'2 '^c) 92hrhc {dabedcre ‘^farefbce'} 
-\-g 3 mrhc {dabed ere farefbee faeefbre) 

Qb^^r^^c {daredbee ~ fare fbee') 

+gemrmcdabedcre - 'igsUlamb- (17) 

which can be readily inverted. These coefficients are to¬ 
tally defined in terms of ha and the parameters of the 
model. 

Now to the fermion determinant related to the inte¬ 
gration over the fermion fields: we expand it using a 
heat-kernel technique that takes appropriately into ac¬ 
count the quark mass difference^ being chiral covariant 
at each order of the expansion [55l - [^ . 

W[Y] = In I det D| = J yp(t) exp , 

d\De = - d'^ + Y, Y = j 7^(5^ -b iibdki.(t)) 

-b cr^ -b {M, cr}-b -b *75 [cr-b M, <()], (18) 

or 

/ 74 

(19) 

i=0 
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where De stands for the Dirac operator in Euclidean 
space. We consider the expansion up to the third modi¬ 
fied Seeley-DeWitt coefficient bi 

bo = 1, bi = -Y, 

b2 = “2 -1" + ^ds)Y, (20) 

with Aij = Mf — Mj. This order of the expansion takes 
into account the dominant contributions of the quark 
one-loop integrals It {i = 0 , 1 ,...); these are the arith¬ 
metic average values h = + Ji{M^) + Ji{M^)] 

where 

00 

Mm^)= I (21) 

0 

with the Pauli-Villars regularization kernel [1^, [s^ 

p{tls3) = 1 — (1 -I- th^) exp(—tA^). ( 22 ) 

which is equivalent to the sharp 4D cutoff regularization 
for the scalar integrals considered. Both terms propor¬ 
tional to hi and &2 have contributions to the gap equa¬ 
tions and meson masses, but only 62 contributes to the 
kinetic and meson interaction terms. By excluding the cr 
tadpole from the total Lagrangian, one obtains the gap 
equations 

h, + [ 3/0 - (3Mf - M") h] = 0. (23) 

where + Mj -|- M^. We now see that (I23p must 

be solved self-consistently with the SPA equations (flfit . 


2 (M„ - + M„M, - Mj] {(j)l + cj)l) 

2 (Md - Msf + MdMs - {(t>l + <(>2)} 




^ah 


(24) 


The kinetic term requires a redefinition of the meson 
fields, 




= go- 


R 
a ’ 


= g(l)a, 5^ 


47r2 


(25) 


to obtain the standard factor 1/4. The flavor and charged 
pseudoscalar fields are related through 



/ 

^/2 

7T^ 

K+\ 


7r“ 

4>d 

V2 

K° 

(26) 

\K- 


<P3 1 

02 / 



and similarly for the scalar fields. In the following we 
concentrate on the diagonal components of (IMl) . which 
according to eq. © induce mixing between the 0,3,8 
field components of (I24|) . in general. Indicating also the 
result of the transformations discussed below in (12^ . (130]) . 
we arrive at the following useful relations among fields 


, \/2(/o + <(>8 ^ , 

03 H-7=- — 03 + b" 


0d = 


-03 + 
2 


73 

7200 


73 


= -03 + Pn 


3^,0-^ = 


(27) 


III. Mixing in the TT'^,r],r]' system 
III a. Choices of representations 

Finally one is ready to combine the terms of the total 
Lagrangian L that contribute to the kinetic terms Lkin 
and meson masses Lmass 


Akin 4“ Yjjiass 

- ^ { [2 (M„ + Mdf - MuMd - M/] {al + aj) 
+ [2 (M„ + M, f - MuMs - Mj] [al + al) 

-\- \2{Md + Ms)^ — MdMg — ((Tg -|- (Ty) 

+ ^ [ol (8M2 - Mj - M/) + (8Mj - - M/) 

+ al {ml - Ml - Ml)] 

+ i [02 {2Ml - Ml - Ml) + 02 {2Ml - Ml - Ml) 
+ 02 {2Ml - Ml - Ml)] 

+ \2{Mu- Mdf + MuMd- Ml 


The neutral physical states 7 r°, 77 , 77 ' are related to the 3x3 
symmetric pseudoscalar meson mass matrix of elements 
Bij emerging in the i,j = {0, 3, 8 } channels of Lmass by a 
sequence of two transformations S = lAV that diagonalize 
it 



^ B33 Bo3 B33 ^ 


03 \ 

( 03 , 00 , 08 ) 5 ' ^5 

Boo Boo Bqs 
y B3S Bos Bss y 

5-15 

00 

Ws/ 


(28) 


first a rotation to the strange-nonstrange basis through 
the orthogonal involutory matrix V 


with 


^ 03 

Vns 
\ Vs 


= V 




(Vs 

0 

V 0 


0 

72 

1 



(29) 


(0? + 0^) 


( 30 ) 
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and then through the unitary transformation U to the 
physical states 


1 

\ rj' 


= U{€i,e2,tp) 


^ </>3 
'Hns 
\ Vs 


(31) 


where 


U = 


/ 1 

— £2 — £1 COS^ 

y —eisin^/> 


£1 + £2 cos ip 
cos Ip 
simp 



(32) 


having a single mixing angle involved in the determina¬ 
tion of the decay constants associated with the r] and rj' 
mesons. We show in the next subsection that our model 
fulfills this condition at the approximation considered. 

At this point one should remind the reader how the 
model’s mri' contains a term related with the topological 
vacuum susceptibility. The generalized NJL Lagrangian 
which combines the 17^(1) breaking by the’t Hooft {2Nf) 
determinantal Lagrangian with the 4q and 8 q interac¬ 
tions has been shown in to be in correspondence with 
the Witten-Veneziano formula which relates , in 
the large Nc limit of QCD with massive quarks, to the 
topological susceptibility xym of pure Yang-Mills, 


The conventional definitions e = e 2 +ei cos ip, e' = £i simp 
for the mixing angles are used in the tables. The unitary 
matrix W has been linearized in the tt *^—77 and Tr^ — rj' mix¬ 
ing angles £ 1 , £2 ^ 0(d), d <C 1. This can be done because 
(p 3 couples weakly to the rjns and states, decoupling 
in the isospin limit, while the mixing for the rj — rj' sys¬ 
tem is strong. Nevertheless we have tested numerically 
the linearization by obtaining also the exact Euler angles 
associated with the transformation (l3T]) . the differ¬ 
ences lying within the one to two percent level for the 
cases studied. 

In the isospin limit U in (1321) leads to the 2x2 or¬ 
thogonal tranformation 

(33) 

with 

(34) 

\ Simp cos Ip J 

We remind that the angle ip is related to the mixing angle 



of the physical states rj, ij' in the singlet-octet renormal¬ 
ized basis states (p^j^, cp^ as ip = 9p + arctan-\/ 2 , with the 
principal value of the angle 9p comprised in the interval 
— (ttM) < 9 „ < ( 7 r/ 4 ), for details please see Appendix B 
of [dll and 

We emphasize that one can freely choose among the 
different orthogonal bases to address the mixing of 
states. The reason it is convenient to adopt the strange- 
nonstrange basis is that it allows to infer whether the 
mixing-parameters are determined in a process indepen¬ 
dent way; in the context of ChPT it has been shown to 
be so if certain OZI violating processes are suppressed, 
[13; [13 (13) and with the exception of those originating 
from topolgical effects due to the U{1)a anomaly. At this 
level of accuracy the decay constants follow the pattern of 
particle state mixing in this basis which is tantamount to 


nip, +mp- 2m\ = --^xym- (3d) 

J TT 

In our model the topological susceptibility is obtained in 
the large Nc limit as the following combination of model 
parameters 

XFM = ^(^)- (37) 

The cutoff A, which is the approximate scale at which 
dynamical chiral symmetry breaking sets in, does not 
appear explicitely in the relation for xym, only hidden 
in the constituent quark mass M. This expression shows 
a judicious interplay of the subleading in Nc counting 
CAi(l) breaking parameter k and the LO in Nc 4q cou¬ 
pling G, which combine in a relation that survives in the 
large Nc limit. This shows that the roles of chiral sym¬ 
metry breaking and the breaking through the Adler Bell 
Jackiw anomaly are intertwined to equip rj' with its large 
mass, which does not vanish in the chiral limit. 

Ill b. Decay parameters 

At the order of the heat kernel considered, we obtain 
the model’s axial-vector current as [13 [13 

71“ = 

itr[({CT^ -k Mg-^,dp,(p^} - (p^})Xa] + 0 (^ 3 ) 

(38) 

in terms of the anti-commutators involving the bosonized 
and renormalized fields cr^, cp^ (| 2 ^ and the constituent 
quark mass matrix M = MaXa- From here it is straight¬ 
forward to calculate the matrix elements in the singlet- 
octet basis 


<o\A;m4>’k>=^r%- ( 39 ) 

One has 

noo _ + Md + Ms _ 22 _ p 33 _ Mu + Md 

^ ~ % ’ ~ “ 2 I 

rss _ A + 4Ms Qg ^ Mu + Md — 2Ms 

65 ’ SV2g 









/03 = y/2f^ 


Mu-Md 

Vdg 

r44 _ r55 _ + Ms g0 _ 77 _ Md + Mg 

2g ^ J - J - 2g 


(40) 

(41) 


In particular one obtains < 0\Ajj^^‘^{0)\n{p) >= iy/2f^p^ 
and < 0 |^ 4 +. 5 (o)|i^(p) >= with U = 

and Jk = ^" 23 ^° ’ order of the heat kernel expan¬ 

sion considered. 

The neutral axial vector currents can alternatively be 
taken in the strange non-strange basis 



for which one obtains the decay constants 

< 0 |yi;^( 0 )|^^ >= {cr,T} = 3,ns,s (43) 


fus 


,3 Mu + Md „„ M., 


25 


^ r = 


j3,ns _ _ jns,s _ q 

2ff 


(44) 


where are short-hand notations for 

f3,3jns,nsjs,s^ The elements of gS]) are collected 
in the following matrix 


A = 



z/tv 

U 


\ 0 0 



(45) 


with 


C. = 


NcIq 2 2 ^ R _ 2 2 ^ 

27 r 2 ^ 3 ^*’ * “ 2^2 ^ 3 ^*’ 






if) ^ai^bj 1 


2 

Cij — ^aiS-bj, (48) 


where are symmetric quantities. The quantities Q 

are 


c. = 2Mf - Mf - Ml iljlk. (49) 

Inserting them in Bi and using the gap equations (1231) . 
one obtains 



(50) 


Similar expressions can be obtained for the but 

they are not needed in the following. The divergence of 
the axial current (1551) reads in this basis 


df,Al = {a^R H- -)d‘^(j)iR - {d'^aM)(j)iR. (51) 

g 

Using now the equations of motion for the atR, (j)iR fields 

0 — d (JiR CidiR “b 

0 = (I 4^iR Bi(j)iR {Cij4*jR “b Cij4'iR') (^2) 

one obtains for the matrix elements 

< 0\df,A^lcl,jR >= + 2C„), (53) 

which encode all the chiral and tb 4 (l) symmetry breaking 
terms. The off-diagonal elements of the inverse matrix 
pertinent to C,ij are calculated to be 


_/V/^ r 3 , 7 is 

with z = = f — marking the departure from 

isospin symmetry. It is of the order of the ratio involving 
^ 0{S), with ff = Mf/g the decay constants 
Tf = diag{fu,fd,fs}- 

According to the idea behind the strange-nonstrange 
basis, the transformation to obtain the physical decay 
constants 

Tp=UJ, P = {7r°,r),r7'}, (46) 

is the same that transforms the states, eg. (1511) . In the 
following we discuss how the obeservables obtained from 
our model Lagrangian fulfill this condition. In order to do 
so, it is convenient to express the meson mass Lagrangian 
Cn of the neutral states in the flavor basis 

i—u,d,s ^ 

+ ^ {^ijO'iRdjR + Qj(j)iR(t>jR) ? (47) 

j—u,d,s I 


Cil = ^{hkK + 2mkK2 + 2gsm^mj), j ^k. (54) 

These contributions violate the OZI rule, the ones pro¬ 
portional to K, k ,2 have their origin in the Ua{^) anomaly. 
The term gg is a contribution of the order of the square 
of the current quark masses. 

For comparison the transition elements of the diver¬ 
gence of the axial current in the strange non-strange basis 
are 

Q) =< 0\df,A1\(l)jR >, i,j = {3,ns,s} (55) 

Ql = ^ib+ + z{b--2(:^d)) 

Qns = ^(.b + z{b'^ + 2(^d)) 

Qr = Y^b-+z{b+-2C^d)) 
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Qnt — + 2(„^)) 


”*" ”*” — Cds)) 


L 


fs 


Q* = - Cds), QL = ^iCus + Cds) 


V2 

Ql = /•^("^ + Css) 




(56) 


with 


— {-^ + C«u)±(-^ + Cdd) 


(57) 


Note that the elements Q®, i,j = {3,ns,s} are not sym¬ 
metric under exchange of {i,j} differing by terms ^ z and 
y — ^ ■ Obviously in the isospin limit the elements 
Qls^QT vanish. In this limit the 
y being a measure for flavor breaking in the light vs. 
strange sector. 

Finally using the relations (|5T|) . (I5^ and (1551) we cal¬ 
culate the vacuum to physical particle transitions of the 
divergence of the axial current, < >, b = 

{3,ns, s}, P = {7r°,y, y'}, discarding terms ^ <5^, and 
are able to show that it fulfills the following relation [2^ 

< Q\d^Al\P, >= (58) 

with T given by (l45l) and M.aa' the physical meson mass 
matrix, which on the requirement of being diagonal de¬ 
termines the mixing angles ip,e,e'. Thus the decay con¬ 
stants transform as the states within our model 
calculations, at the order of the heat kernel considered. 
Higher order terms involve derivative interactions, which 
are likely to change this behavior. 

We obtain the foilwing relations at 0{6) 


E £3 £3 £2 \ £ns £ns £2 

JPJP — / tt ’ 2^ JP JP — /tt 

-7r^,r),r)' P=7 t^ ,r},r}' 

E fpfp = fs^ 

- 7 T^ , 77 , 77 ' 

(59) 


for the diagonal elements, and 


case. In the isospin limit the physical states P = rj^rj' are 
now only mixtures of the nonstrange and strange compo¬ 
nents, as follows from (l33)) . This case has been considered 
in [62j| (where a detailed discussion of the mixing scheme 
in connection with the singlet octet basis in the isospin 
limit is also presented), 

(OI-^m(O)I-P(p)) = i/pP/., (i = ns,s). ( 61 ) 

The couplings can be represented in a way which is sim¬ 
ilar to the one of Leutwyler - Kaiser , who intro¬ 

duce two mixing angles •dns , '&s 

\f\=i \ ( frisCOSdns -fsSm^s \ 

^ v/v j yfnsSin-dns fs cos-ds j 

(62) 

Our calculations show that within our model, with M = 
Mu = Md, 


= — cos f), /® =-^ sin V', 

5 ' 9 

/”/ = — sin fC,' =— cos Ip. (63) 

9 9 

One thus obtains the relation for the mixing angles 

E fPfp = rrsin[dns-^s)=0. (64) 

P—7).,r)' 


This follows as the basic parameters fns, fs, '^ns,'&s of the 
matrix {fp}, expressed in terms of model parameters (in 
the approximation considered) 


fns = - = U, fs= - 1p='dns='&s- (65) 

9 9 

There is a direct relation between the common mixing an¬ 
gle -dns = '&S and the OZI-rule which has been discussed 
in (63,(55. The model predictions agree well with the 
general requirements of chiral symmetry following from 
chiral perturbation theory (ChPT), although the results 
differ already at lowest order. For instance, we have 


f! = Pfl -fl + 


(M, - Mf 

2y2 


( 66 ) 


E /pV|.= E /p/p = o> 

P—TT^ ,r},r}' P— 7 T^,r},r}' 

E fpfp = ( 60 ) 

P^TT^, 77 , 77 ' 

for the crossed terms. The last relation is a consequence 
of (HSl) . The vanishing of the other crossed terms indi¬ 
cates that just one mixing angle is present for the y — rj' 
mixing in the strange nonstrange basis at this order of 
approximation. This can be seen in a simple way in the 
isospin limit, and the argument is the same in the general 


where barred quantities are a reminder that they are 
taken in the isospin limit. The first two terms of this 
formula are a well known low-energy relation which is 
valid in standard ChPT. In the y'-extended version of 
ChPT there is a OZI-rule violating term in the effective 
Lagrangian, which contributes as pAi to r.h.s. of (l66ll . 
We have instead the term {Ms — M)^/(2y^). Of course, 
in our case the origin of this contribution is related with 
the SU{3) flavour symmetry breaking effect and does not 
have impact on the deviation from a single mixing angle 
in the strange-nonstrange basis reported as consequence 
of the OZI-rule violating parameter Ai. 
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Away from the isospin limit one obtains 

= + (67) 

with a correction of order d in the term, as compared 
to (l66)) . The numerical values that we obtain are z ^ 
4 X 10“^, which amounts to a small correction of ^ 0.8% 
in the term. 


For completeness we indicate the decay constants and 
mixing angles in the singlet-octet basis in the isospin 
limit, used in Table dm, for details see please [62|, 


^ M 

^0 = Ip — arctan{\ 2^^), 

9s = tp — arctan{\/2^^) 
Ip = 9p + arctan'/2 


2 _ 2/i + /^ p M, 

° _ 3 _ 6 ^ M 

.2 _ 4/i-/2 , {Ms-Mf 

“ 3 + 3g2 


IV. Results and further discussion 


( 68 ) 


(69) 


The model has 15 parameters, 4 couplings G,K,gi,g 2 
associated with L^t in 7 non-vanishing couplings 
K 2 ,g 3 ....gs in in the ESB sector ([4]), the cutoff A, and 
the 3 current quark masses. Before running a fit it is 
convenient to understand to which parameters the dif¬ 
ference of the light quark masses {rriu — rrid) is most 
sensitive. For that it is instructive to look at the ma¬ 
trix components dm in the strange-nonstrange basis, 
which vanish in the isospin limit. As mentioned before 
these elements belong to the SPA contribution to the 
meson mass Lagrangian, see last line in (I24L which is 
the part that carries the full information on the model 
couplings. Note that the heat kernel contribution to the 
meson masses, represented except for the kinetic terms 
in the remaining of expression (1241) , only depends on the 
cutoff A of the li quark integrals. The dependence on 
the model couplings of the heat kernel contribution only 
enters implicitly through the constituent quark masses, 
via the gap equations (ESI) which are solved selfconsis- 
tently with the lowest order SPA equations ((T51) . Defining 
toa = \{md-mu), rriY. = \{md + mu), = ^(hd-K) 
and hs = \{hd + hu), one has for the inverse matrix 
elements of the (pij matrix, i,j = {3,ns, s} 

(Cs.ns)^”^^ = ^[/iA(2g2^S + gs^s) 

+mA{g3h^ - 2(^5 - 36 + 238 ) ws )] 


(C3..)(-') 


1 

271 


{Hak + 2mA{K2 


gsrris)) 


(c 

\Sns,s ) 


1 

2 ^ 


(/lEK -I-2(6)2 + 38 TOs)toe). (70) 


The first two elements vanish in the isospin symmetric 
case. The third element relates to the mixing in g — g' 
and is non-zero in this limit. In this basis one sees that 
the mixing in the 3, ns sector involves other couplings as 
in the strangeness related sectors, in contrast to the off- 
diagonal matrix elements in the u, d, s basis, ([S3]). This 
sets new constraints on these couplings, as compared to 
the isospin symmetric case. It follows that the interplay 
of these parameters (which is conditioned by the fits), not 
the actual magnitude of each of the terms, is relevant to 
obtain the size of isospin corrections. We remark that in 
isospin breaking effects are absent, and that in 
the diagonal elements (not shown) they are overshadowed 
by the presence of mg, ms, fig, /is. 

The current quark mass dependence in these expres¬ 
sions enters together with the ESB couplings. In their 
absence the effects of ESB come only through the differ¬ 
ence in the light condensates /ia which do not vanish if 
the conventional QCD mass term has ruu ^ md values, 
and if the couplings k and 32 are not zero (if they also 
vanish, only the heat kernel contribution to the meson 
mass matrix carries the effects of isospin breaking). The 
coupling K is strongly correlated with the g — g' mass 
splitting and it enters in the corresponding (Cns.s)*'”^^ 
matrix element as factor of /is which remains approxi¬ 
mately constant. Thus this parameter is not expected to 
vary much in the fit of isospin breaking effects, which has 
been also verified numerically. 


Before showing the results for isospin breaking, we dis¬ 
play relevant model observables in the isospin limit, in 
comparison to other approaches. We consider the cases 
in which the ESB terms are present in the interaction 
Lagrangian, sets (a,b) in the Tables Ill lIVl and compare 
with the parameter set (c) in which explicit symmetry 
breaking occurs only through the LO current quark mass 
term. Table I indicates the mass spectra of the low ly¬ 
ing pseudoscalar (and scalar meson nonets, see caption, 
for sets (a,b)) used in the fit of parameters. Sets (a,b) 
have different mixing angles for the pseudoscalars as well 
as for the scalars, leaving however the mass spectra and 
weak decay constants inchanged. This is possible 

because the fit is done simultaneously in both sectors, 
leading to a readjustment of parameters. A comment 
on the scalar mass spectrum: we were able to obtain a 
reasonable simultaneous fit for the a mass and its large 
decay width, at tree level of the mesonic Lagrangian. It is 
our understanding that this results partly from the fact 
that already at mesonic tree level there are signatures 
of quark-antiquark as well as of admixtures of 2 quark- 
2 antiquark states present, from the underlying multi¬ 
quark Lagrangian, which contribute in the long distance 
asymptotics to the internal structure of the mesons as 
well as to the mesonic interaction terms. In many other 
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approaches the more complex quark structures enter only 
through the explicit inclusion of meson loops, tetraquark 
configurations, and so on 

Table II shows the current and constituent quark 
masses in the isospin limit and the model parameters 
which do not break explicitly the chiral symmetry, the 4- 
quark coupling G, the 6 -quark t Hooft determinant cou¬ 
pling K, and two 8 -quark couplings gi,g 2 , of which gi 
is OZI-violating. In set (c) we put gi = 2500 without 
loss of generality, since a change in gi can be counter¬ 
balanced by a change in G leaving all other parameters 
and observables unch ang ed, except for the low lying a- 
meson mass, see e.g. [^, [^, which does not affect the 
pseudoscalar characteristics considered here. Table III 
displays the couplings related with the ESB interactions. 

In Table IIVI further properties of the rj — rj' system are 
presented, the angle dp in the singlet-octet basis and the 
related angles 9o,9s, as well as the weak decays /o, /s) see 
eqs. (IB51) . (I5^ obtained following the methods in sections 
3.2 and 3.3 ofjo, and compared with the results of other 
approaches [22j| and references therein, [t^ , , [131 , and 
of Fade approximant method to pseudoscalar transition 
form factors [i3|, [Hi, where in the Table ITVl ( 77 ) and ( 77 ') 
refer to the asymptotes taken for the fits, and the results 
obtained from BABAR data [s^l by the same authors 
(see discussion in [HI concerning the BABAR ’’puzzle”). 

This comparison shows that set (a) yields results quite 
close to the chiral perturbation analysis of [zl] and fol¬ 
lows the general trend of the other references in the table, 
with exce ptio n of the sum rules approach of [131 , and the 
fits from |83| . Our values for 9o,9s are smaller, respec¬ 
tively larger than the ones presented in 22 |, which is 
probably due to the different way in which the U{1)a 
anomaly is treated. One further notices that an increase 
in 9p leads mainly to an increase in the 9o angle, com¬ 
paring sets (a) and (b), which have similar theoretical 
input. 


Turning now to isospin breaking, in the numerical fit 
we keep the cutoff A, which sets approximately the scale 
of chiral symmetry breaking, close to the isospin limit 
value, and take for the strange current quark mass the 
value ms = 95MeV. We do not consider the scalar spec¬ 
trum in this case, the related isospin breaking effects will 
be addressed in a future work. Then the general case, 
with the ESB terms, involves solving self-consistently a 
system of 13 equations, the three gap equations (IMl) sub¬ 
ject to the three SPA conditions ®,6 equations for the 
meson mass matrix elements, that is 3 for m^o, m^, 
and the remaining 3 for diagonalization, one equation for 
the kaon mass, 2 equations for the weak decay constants 
Jtt, fK, one equation which fixes the ratio r = of 

current quark masses. We vary externally the values of 
the mixing angles and search for the parameters 


that lead to the best fit of the mu,md current masses. 
The result is indicated in Table VII, using empirical in¬ 
put and couplings shown in Tables V and VI (the empir¬ 
ical splitting of the charged multiplets cannot be repro¬ 
duced, since electromagnetic effects are not taken into 
account). In presence of ESB (sets A and B) with the 
ratio r kept equal and close to its empirical value, the 
ratio ^ is well reproduced in comparison with the liter¬ 
ature shown in Table VIII. The main observation is that 
this ratio is reduced by ^ 40% compared to the model 
variant without ESB interactions, set C. In the latter 
case one does not obtain 4 nor TO„,md close to the em¬ 
pirical values. Neither this ratio nor the light current 
quark masses get improved in set C by reducing r down 
to 0 . 2 , the main consequence being a drastic change in 
the values for the mixing angles, which get reduced to 
e = 0.0119, e' = 0.00135. We note that it is not possi¬ 
ble to get a better fit for m,;, mK and fx in absence of 
ESB interactions as the one shown in set C in Table V. 
Contrary to this, the individual values for e and e' for 
sets A,B are in good agreement with the ones indicated 
in 1 ^ . 25 1 and the corresponding current quark masses 


are very close to the quoted values m„ = 2.15(15) MeV, 
TOd = 4.70(20) MeV, [s^. 

Regarding table VIII that collects values obtained 
in the literature, within different phenomenological ap¬ 
proaches, as well as in experiments, a comparison of the 
different values has to be done with care, for a careful 
and detailed discussion see [131. We note in particular 
that in the experimental value [83| the mixing 7 r° — rj' 
has not been taken into account and that in the ChPT 
result [ 13 ] the 77 ' is considered as a background field. 


V. Conluding remarks 


We conclude that the explicit symmetry breaking in¬ 
teractions of the generalized NJL Lagrangian considered 
are crucial to obtain the phenomenological quoted value 
for the ratio We obtain values for the e mixing angle 
which lie within the results discussed in the literature. 
Unfortunately the value for e' is much less discussed. We 
obtain e and e' reasonably close to the ones indicated in 
[ 131 , [ 13 | for current quark mass values in excellent agree¬ 
ment with the presently quoted average values. The cor¬ 
responding sets (A,B) are the ones which also yield the 
best fits to other empirical data within the model vari¬ 
ants. 
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